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ABSTRACT 



Context. We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel 
Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT 
Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular 
fs^ , cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. 
pTT ' Aims. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to 

. ' give an inventory of the gas lines present in its far-infrared spectrum. 
: <""| Methods. The 69yum feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 
O i, 32 emission lines from five gaseous species and measured their line fluxes. 

I ■ Results. The 69yum emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises 

from ~200K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have 

.£3 [ few defects and contain at most a few percent iron by mass. Forbidden line emission from [C n] at 157 fim and [O i] at 63 and 145 pm, most likely 

C/3 due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational 

transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J = 22 - 21 and T ~800 K for higher transitions. 
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f--» 1. Introduction ratio. iBouwman et alJ (l2003h suggested an inner c avity in the 

tt|- _. „ , , , . , , . disk. Later observations confirmed this gap (e.g., Gradv et al. 

^ , Circumstellar disks around young stars are the birthplaces of [ 2QQ5 . Benisty -J^i] [20H). This remarkable feature and the 

. ■ planetary systems. To understand planet formation, it is vital lon it of me HD 100546 disk int to a lanet 

jn , to study the processes that govern the evolution of gas and fa ^ inner 1Q ^ maki k a ime for detailed smdies 

O" dust in these disks. PACS provides unique information in this „, „ TTr . * nn c,<- , • , ~ , , on 
,. , , , , , . j, j. . j, , , T „, ,. , „ The system of HD 100546 has a rich 2.4 - 180 um spec- 
field through observations of far-infrared (IR) solid-state tea- / , . ., T „ „ k „ , c .. . , | lr ,r, k i_ • ;•, 
t— I . , . ° . , , .. v ' ^ , . Uum as observed with ISO (Malfaitetal. 1998), showing a stak- 
. . tures which are particularly sensitive to temperature and ele- .... . . , .Jr TT , _ — 17^ 7^ — ri 

>. . . .„ r,»^n • 11 • j j 1 ing similarity with that ot the comet Hale-Bopp (Crovisieret al. 

mental composition. Moreover, PACS is well suited to study the 1 ,^,-4 , i T ,, — — n 1^77771 . .,,, „ ., c 

, , ,„„,. . . . J , , 1997). Juhasz et al. (2010) investigated the Spitzer spectrum ot 

k> warm gas of a tew 100 K, complementary to the hot gas probed ~r~ — 7T~ . — . . , , . r , ., TO „ 

\^C ■ . ° TT , , , , , , , .„. , , the source. Strong forsterite emission is observed in both ISO 

r.^ in near-IR and the cold gas observed at millimeter wavelengths. , c .. , . . t . ». . ,. . , , „. , . . 

$-h »••,,., • • , tt , • and Spitzer data. A continuum fit indicated crystalline dust at 

^ • on iS mtenS1V ^ f A U A d f/ P r '"T n " SeqUen f StarlS g 210K and 50K (Malfait et al 1fl998h . 

B9.5Vne star HD 100546. While signs ot ongoing accretion _ Tr ,_. : ,. . r ^, . ,„ , 

. ,. , ., , , . j. l r. , ., " i i Unn/il T tt, »■ « a The ISO spectrum shows a strong line ot [O 1] 63 um and a 

indicate the stars infancy (Deleuil et al. 2004), the estimated . r _ n \ -„ ,. T „ TT > \ n -^ ■ — 7^1 1^™^ 

c lnl >T , .. 7. TTV 1 -n. a- 1 weaker [C 11] 158 um line. In the near-IR, Brittain et al. 2009) 

age of lOMyr makes it unusually old for a star with a disk .1 "-: — ; T ^ t: — - 1 rrr— -1 ' : — ' *—. — " 

IT — a ai ^ l ll inmK rminnfK' u /1m \ a an d van der Plas et al. (2009) observed CO ro-vibrational lines 

(vanden Anckeret al. 1997). HD 100546 is nearby (103 pc) and .' — — ' , — — " . ,„„„„ ,„ , „„ 

.. ,, Tri „ ai- t.4. ■ u 1 j i.u c with a rotational temperature of ~ 1000 K and found no CO gas 

optical/neai--IR scattered light imaging has revealed a wealth of . . . , , A r _. . , . . r _. . , _ , 

f^ . ., ,. , , ° rp, — . P i l hr>r>i L 17 TTI in the inner dust cavity. Observations of optical [Oil and Balmer 

structures in the disk (e.g., Gradv et al. 2001; Augereau et al. ,. . . , , , , , , , ■ 

hnni L U j-i — mfinmi d a ^u a * tt> line emission, on the other hand, demonstrate th at the inner dis k 

2001; Ardilaet al. 2007). Based on the mid- to near-IR excess . ' , , . , , Vr— ; — 2 1 : — \ — TH ^™A . 

is not completely devoid ot gas (Acke & van den Ancker 2006). 
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that was obtained within the Herschel key program 'Dust, Ice 
and Gas in Time' (DIGIT). 



2. Observation and data reduction 

We used the spectral energy distribution (S EP) mode, cover ing 
the wavelength range from 55 to 210 urn (see |Poglitschl2010l for 
details of the instrument and calibration). Four up/down scans of 
the grating for two nod positions were performed with a chopper 
throw of 1'. Our data were processed in HIPE 3.0, provided by 
the Herschel science center, up to level 1 . We did not apply the 
integration ramp linearity, cross-talk and response drift correc- 
tions as these were not yet scientifically validated. 

After level 1 processing, we extracted the spectra for the 
two nod positions separately from the central spaxel (spa- 
tial element). The wavelength-dependent point spread function 
(PSF) correction was taken from the PACS calibration document 
PICC-KL-TN-04 1 . We first uniformly rebinned the spectra, us- 
ing an over-sampling factor of 1 relative to the spectral resolu- 
tion of the instrument and an up-sampling factor of 2, before 
combining the spectra of the two nod positions. The current flux 
calibration and relative spectral response function (RSRF) were 
derived from ground-based measurements. The correction fac- 
tors, determined from the first in-space calibrations were applied 
(for details see the PACS calibration document PICC-KL-TN- 
04 lQ). We estimate that the absolute flux calibration is accurate 
within 30%. The RSRF is valid for spectral features not much 
wider than about 1 fim. For broader features, the RSRF is not 
yet well calibrated. Therefore, we will focus on spectral fea- 
tures narrower than about 1 fim. The wavelength calibration is 
accurate down to the spectral resolution as can be seen from the 
gas-phase lines, most of which are found within 0.02 fim of their 
laboratory wavelength. 



3. Results and Discussion 

The PACS data confirm the tentative ISO detection of the 69 fim 
forsterite band first reported by Malfait et alj (11998) and deliver 
a band profile of very high quality (see Fig. [TJ. The feature is 
characterized by a signal-to-noise ratio of about 10. In addition, 
we find a wealth of narrow molecular gas lines in the spectrum 
of HD 100546, coming from CO, H z O, and OH, and [C n] and 
[O i] fine-structure lines (see TableQ]). 

In this section we present an analysis of the 69 fim forsterite 
peak in terms of dust temperature and composition. In addition, 
we provide a compilation of the detected gas lines in the far-IR 
spectrum of HD 100546 and derive the CO-rotational tempera- 
ture. A full spectral decomposition, including a search for shal- 
low and wide features (e.g. ice and hydrated silicate features), 
will be presented in a follow-up paper. 



3.1. Forsterite 

The '69 fim' band is a characteristic signature of f orsterite - 
the Mg-rich end member of olivines dHenn ing 2010). In order 
to determine the position and width of this feature, we fitted a 
Lore ntz profile and a lin ear function for the continuum (compare 
e.g.. iKoike et al.l 120 06) in the wavelength interval 68-71 fim. 
The emission peaks at 69.198 + 0.008 fim and has a FWHM of 
0.64 ± 0.03 fim. 
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1 http://herschel.esac.esa.int/Docs/Calibration 
/PACS_SpectroscopyPerformanceCalibration_ HMar2010.pdf 



Fig. 1. The continuum-subtracted forsterite emission at 69 /im (his- 
togram) overplotted with the best-fit model (red line) which is given 
by eq. (TJ. The contributions of the different temperatures are indicated. 
The dotted vertical line locates the peak position of 70 K pure forsterite 
and the dashed one the position in the PACS spectrum, fitted with a 
Lorentz profile, which is not shown here. 



The 69 fim band pr ofile depends on bo t h the t emperature and 
the iron content (e.g., IKoike et al.l l2006l 120031) of the olivine 
crystals. Increasing the dust temperature leads to a shift of the 
peak position towards longer wavelengths and to an increase of 
the FWHM (Fig. |2). Increasing the iron content in the olivine 
crystals has a very pronounced effect in shifting the band to 
longer wavelengths (Fig. |3). Finally, the FWHM will increase 
with the number of lattice distortions (Koike et al. 2 01QI) . 

Based on the peak position of the band, we can immediately 
conclude that the olivines cannot contain more than 3>-A% iron as 
this would cause the peak to appear at significantly longer wave- 
lengths than observed. The PACS observations provide much 
better constrai nts on the iron con tent than the shorter- wavelength 
Spitzer data (Juha sz et al.ll20101) because silicate bands in that 
spectral range are less sensitive to a change in Mg/Fe ratio. 
Furthermore, the small FWHM of the ba nd profile indicates 
crystals with few defects (IKoike et al.l l2010). The small FWHM 
also shows that the dust emission comes from a narrow range of 
temperatures. 

Determining the temperature of the dust that produces the 69 
fim feature from intrinsic properties of this band is more difficult 
because even a small amount of iron could balance a lower tem- 
perature. Assuming pure forsterite results in a dust temperature 
of about 200 K (Fig. 0. To constrain the amount of colder dust, 
we fitted the observed profile with a weighted sum of laboratory 
data (Sutoetal. 2006) at different temperatures, ranging from 
50-300 K (Fig.[TJ. The measured flux F(,^ 71flm is given by 

5 

Fet-nprnW = F cont (A) + J] w t • K(r dm , A) ■ B A (A, T' du J. (1) 

Here F cont is the local continuum, and the k values are the mass 
absorption coefficients, computed from the optical constants of 
the Suto et al.l ([2006) forsterite sample, at the different dust tem- 
peratures T\ , Finally, the w,- values are the relative weights 
of the flux at the different temperatures. In addition, we had to 
specify a dust shape distribution in order to calculate the mass 
absorption values. We used the DHS model (distribution of hol- 
low spheres, Min et al. 2003) with a grain size of 0.1 fim, which 
has been foun d to be a good rep resentation of observed silicate 
:alJ201(l) . 



profiles (e.g. Jjuhasz et i 
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Fig. 2. FWHM of the forsterite '69yum band' plotted over the peak po- 
sition, fitted with a Lorentz profile. The square is our HD 100546 data 
(size indicates error). For comparison, data from Koike et al. (2006, 
asterisks: 'Kyoto sample', triangles: 'Jena sample'). ISuto et all (2006. 
crosses, DHS with 0. 1 pm. grainsize) are shown. 
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Fig. 3. Effect of iron-rich olivines on the peak wavelength of the 
'69yum' band at room temperature. T he offset in w avelength is rela- 
tive to pure forsterite. The data are from lKoike et alJ ( 120 03) and the line 
is the best-fit linear function. 



Under the assumption of pure forsterite, no more than a few 
percent of the dust emission can be coming from the 50 K com- 
ponent. This could be up to 9%, if grains of ~5 pm radius are 
assumed. Such large grains can only exist in the deeper layers 
of the disk, where they are invisible to Spitzer, as lJuhasz et alJ 
(1201 Ol) rule out grains > 1 pm. The bulk of the emission, under 
the assumption of pure forsterite, comes from ~1 50-200 K. 

However, a temperature of 150 to 200 K for the 69/im band 
would lead to tremendously higher fluxes in the 20^1-0 pm re- 
gion than observed by ISO and Spitzer. Additionally, the pre- 
dicted relative strength between the 24- and 33 pm forsterite 
bands would be inverted compared to the observations. In fact, 
temperatures of onl y 70K (Bouwman et al. 2003) or even 50K 
(Malf aitetalJI 19981) were derived from the ISO data. 

These results can be reconciled in two ways. First, the dust 
producing the 69 pm feature is at 70-50 K. The shift from the 
band position in the 70-50 K laboratory data to the one observed 
by PACS could be caused by an admixture of iron. Given the 
information in Fig. [3] we need at least 2 but no more than 3% for 
such a shift, considering that the data were taken at ~300 K. 

An alternative explanation is based on the optical depth of 
the disk which decreases with increasing wavelength. The r = 1 
surfac e is deep into t he disk at 69 pm, even reaching the mid- 
plane (lMuldersll20ioh . The disk of HD 100546 is a special case 
in this respect: because of its large inner hole, the midplane tem- 
perature at 13 AU is much higher (~150-200K) than in disks 
without a hole so that the far-infrared features can be produced 
in heated layers close to the midplane from which the radiation at 
70 pm can escape. In the most extreme case the emission could 
arise from pure forsterite of 200 K temperature, located in the 
midplane at 13 AU. 

In this scenario the location of the forsterite near the inner 
edge of the disk supports the crystallization by collisional cas- 
cades or shocks induced by the suggested giant planet. 



3.2. Gas lines 

We identified 32 gas lines (Table [TJ from the ISO line listo 
CO lines up to the 7=31-30 transition were detected as well 
as H2O and OH lines. Selected molecular lines are shown in 
Fig-El Furthermore, we confirm the IS O detections of [Oi] and 
[Cn] fine-structure lines (Lorenzetti et al. 2002). The integrated 
line fluxes, obtained from fitted Lorentz profiles, are in Table 
Q] The [O 1] line fluxes agree within 30% with those seen by 




wavelength [micron] 



Fig. 4. Selected gas lines from the PACS spectrum of HD 100546. The 
continuum-subtracted, normalized fluxes are plotted as histograms. The 
continuum has been fitted with a linear function of the wavelength and 
the peaks with Lorentz profiles. The second peak in the H 2 90yum 
panel is CO J =29-28. 



http://www.mpe.mpg.de/ISO/linelists/index.html 



ISO (lLorenzetti et al.ll2002h . while the [C 11] line flux is about a 
factor of 3 weaker. Significant flux outside the pointsource PSF 
is detected for the [Cn] line, which is not seen for any other 
line, suggesting that extended emission dominates the larger ISO 
beam. Through comparison of the flux ratios \Q i\m/\Q 1T 45 anc ' 
[O i]63/[C n] 158 with the models from lKamp et al.l d2010b a gas 
mass of ~ 10~ 3 Mq is derived. 

According to our knowledge, the detection of H2O is one of 
the first in the disk of a Herbig Ae/Be star. The emission seems 
to emerge from optically thick gas. A possible source for the 
H2O could be sublimation of water ice. A detailed analysis will 
be presented in a future paper. 

Our focus here is on the newly detected CO lines. In Fig. 
[5] we show the rotational diagram of the CO lines. Although 
a single-temperature 7 = 580 ± 14 K provides a reasonable fit 
(reduced x 1 - 15) to the data, a better representation could 
be obtained by a two-temperature model (reduced^ 2 = 1.8). 
The transitions up to 7=22-21 are fitted by T - 300 ± 12 K. The 
higher rotational transitions require a hotter gas component of 
800 + 100K providing the connection to the high-temperature 
gas in the inn ermost part of the disk detected by near-in frared 
spectroscopy ( Britt airi et al.ll2009l : lvan der Plas et al.l2 009). This 
inferred temperature range is consistent with CO in a heated sur- 
face layer of the disk as found in disk models that do not assume 
equal gas and dust temperatures but include an explicit treatment 
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Species 




Transition 


Peak 

w 


Lineflux 

10-' 7 [W/m 2 ] 


[Oi] 




3 P, - 3 P 2 


63.19 


554.37 ± 5.0 


[Oi] 




3 Po- 3 P, 


145.54 


35.70 ± 1.3 


[Cn] 




P3/2 _ Pl/2 


157.75 


31.87 ± 1.0 


CO 




7 = 14-13 


186.01 


7.42 ± 0.9 


CO 




7 = 15 - 14 


173.64 


11.50 ± 0.8 


CO 




7=16-15 


162.82 


8.25 ± 0.9 


CO 




7=17-16 


153.28 


10.49 ±0.8 


CO 




7=18-17 


144.79 


9.85 ± 1.0 


CO 




7=19-18 


137.21 


8.86 ± 0.7 


CO 




7 = 20 - 19 


130.38 


7.10 ± 0.7 


CO 




7 = 21-20 


124.18 


5.94 ± 0.6 


CO 




7 = 22-21 


118.59 


5.92 ± 0.9 


CO' 




7 = 23 - 22 


113.47 


4.76 ± 1.2 


CO 




7 = 24 - 23 


108.77 


10.56 ± 1.3 


CO 




7 = 26 - 25 


100.47 


12.71 ± 2.4 


CO 




J = 21 -26 


96.78 


10.20 ± 1.4 


CO 




7 = 28 - 27 


93.35 


10.00 ± 1.0 


CO 




7 = 29 - 28 


90.17 


11.83 ± 1.6 


CO 




7 = 30 - 29 


87.19 


8.46 ± 1.0 


CO 1 




7 = 31-30 


84.40 


4.90 ± 1.5 


OH 


3/2, 


9/2 + - 7/2" 


65.29 


10.00 ± 1.5 


OH 1 


1/2. 


7/2- - 5/2 + 


71.17 


9.17 ± 1.4 


OH 1 


1/2, 


7/2+ - 5/2" 


71.22 


10.10 ± 1.1 


OH 1 


3/2, 


7/2 + - 5/2" 


84.42 


14.96 ± 1.7 


OH 


3/2, 


7/2- - 5/2+ 


84.59 


14.96 ± 1.5 


OH 


1/2, 


5/2 + - 3/2" 


98.73 


8.74 ± 1.5 


OH 


1/2. 


3/2- - 1/2+ 


163.41 


2.60 ±0.6 


P-H20 




322 - 2n 


90.00 


14.32 ± 1.2 


0-H20 1 




4l4 - 3()3 


113.53 


3.64 ± 0.9 


0-H20 




330 _ 321 


136.48 


4.52 ± 0.7 


P-H20 




3 13 - 2()2 


138.54 


3.17 ± 0.7 


o-H20 




2l2 _ loi 


179.54 


5.53 ± 0.8 



1 deblended lines, CO not used in rotational diagram 

Table 1. Peak wavelengths and line fluxes for the identified gas lines 
between 60 and 190yum. Only the formal errors are given, not including 
systematic errors due to the PACS spectral response function. 



of the heating and cooling processes (e.g., Kamp & Dullemond 
2004). The 300 K component could also emerge from the region 
in the midplane close to the gap. 

According to lBrittain et al.l(l2009h the 4.7 yum CO fundamen- 
tal emission shows both a UV fluorescent and a thermal com- 
ponent. The latter was fitted wit h a column density of at least 
2.3 • 10 15 cm 2 (for discussion, see lBrittain et al.l20 09f). The sub- 
millimeter data give a typical CO column density of 10 18 crrr 2 
(Panic et al. 2010). From Fig. [5] assuming that the emission orig- 
inates from the area between 13 and 100 AU, we derive column 
densities of 10 16 to 10 17 cm -2 in both models. These values in- 
dicate that PACS probes the part of the disk between the near-IR 
and submillimeter data. 



4. Conclusions 

The PACS data provide the strongest constraints yet on the 
composition of olivines in a protoplanetary disk. We found the 
olivines to be extremely iron-poor (less than 3-4% iron). Our 
observations can be modeled with dust at ~70 K, consistent with 
ISO results, but requiring at least 2% iron. On the other hand the 
determination of the temperature of the forsterite based on in- 
trinsic features of the 69 //m band instead of using a continuum 
analysis offers a second option: The emission may emerge from 




3000 



EJK] 



Fig. 5. Rotational diagram of the CO lines i n HD 100546 (crosses) , 
overplotted with our models, using eq. 2 from lusttanont et al. (2000). 
A single temperature fit (dashed) gives reasonable results, but a two- 
temperature approach (solid) yields a better description of the data. CO 
7=23-22 and 7=31-30 are blended with H2O and OH respectively and 
were not used in the fit, nor are they shown here. 



pure forsterite at 200 K and 13 AU close to the midplane. This 
region is optically thick for shorter wavelengths. 

The far-infrared spectrum of HD 100546 contains a wealth 
of molecular gas lines including CO, H2O, and OH. We found 
CO transitions up to 7=31-30 which arise from gas in the tem- 
perature range of 300 to 800 K and appear to sample gas in the 
surface layers of the disk at temperatures and surface densities 
intermediate to those probed by submillimeter and near-infrared 
CO measurements. 
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